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ABSTRACT 

We present a new model of molecular cloud cores that originate from filamentary clouds that 
are threaded by helical magnetic fields. Only modest toroidal fields are required to produce 
elongated cores, with intrinsic axis ratios in the range O.f ~ q ~ 0.9 and mean projected axis 
ratios in the range 0.3 ~ (p) ~ 1. Thus many of our models are in good agreement with the 
observed shapes of cores (Myers et al 1991, Ryden 1996), which find axis ratios distributed 
about the mean values (p) 0.5 — 0.6 and q w 0.4 — 0.5. We find that the Bonnor-Ebert critical 
mass is reduced by about 20% by the helical field threading our models. We also find that 
models are generally most elongated when the mass is significantly less than the Bonnor-Ebert 
critical mass for unmagnetized cores. 

Subject headings: ISM: magnetic fields - ISM: clouds - MHD 
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1. Introduction 

There is considerable observational evidence that most molecular cloud cores and Bok globules are 
elongated structures that may be approximately prolate in shape. Myers et al. (1991) measured the 
projected axis ratios of 16 cores in line intensity maps of NH 3 , CS, and C ls O and found that the mean 
projected axis ratio (p), defined as the ratio of the minor to major axis, is in the range 0.5 — 0.6. They 
assumed that the cores are randomly oriented on the sky and calculated the allowed values of the intrinsic 
axis ratio q. Their analysis showed that modest elongation with q = 0.4 — 0.5 can account for the projected 
axis ratios if the cores are intrinsically prolate, but rather extreme flattening, with q < 0.1 — 0.3 is required 
if they are intrinsically oblate. Moreover, 6 of the 16 cores appear to be inside of and aligned with larger 
filamentary structures of lower density observed in CO or dust extinction. They argued that cores must 
almost certainly be prolate in such cases, which are actually quite common. The Schneider and Elmegreen 
(1979) Catalogue of Dark Globular Filaments shows many other examples of elongated cores contained 
within filamentary structures. 

Rydcn (1996) extended the pioneering work of Myers et al. (1991) by analyzing the distribution of 
axis ratios for a large sample of Bok globules and cloud cores taken from several independent data sets (See 
Rydcn for references). She found that randomly oriented oblate cores are inconsistent with the observed 
axis ratios for most samples to a high confidence level (> 99% in some cases). On the other hand, she found 
that prolate cores are almost always consistent with the data. Her analysis predicts that the mean intrinsic 
axis ratio for cores is typically 0.4-0.5, in agreement with Myers et al. (1991), but also that there should be 
significant numbers of cores that are either significantly more elongated or more nearly spherical. 

If cores are truly prolate, then several interesting questions are raised about the non-isotropic forces 
responsible for their shapes. It is well established that rotation is negligible compared to magnetic and 
gravitational forces (Goodman et al. 1993), and tends to make cores oblate in any case (Tomisaka, Ikeuchi, 
and Nakamura 1988a, hereafter TIN). It is known, however, that large scale magnetic fields play a central 
role (Myers and Goodman 1988a, b) in supporting cores against self-gravity. Considerable theoretical 
effort has gone into producing magnetized, non-spherical models of cloud cores (Mouschovias 1976a, b, 
TIN, Tomisaka, Ikeuchi, and Nakamura 1988b, 1989, 1990 (hereafter TIN88b, TIN89, and TIN90), and 
Tomisaka 1991). Most of these models assume that the magnetic field threading the core is purely poloidal. 
However, purely poloidal fields result in oblate structures which, as we have seen, do not seem to explain 
the observations. 
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We show that prolate models of cores, whose shapes are in better agreement with the observations, 
are readily produced by including a toroidal component of the magnetic field, so that the field is helical in 
general []. The basic point is that the radial pinch of the toroidal field helps to squeeze cores radially into 
a prolate shape, while helping to support the gas along the axis of symmetry. Tomisaka (1991) finds that 
toroidal fields can result in moderately prolate equilibria with intrinsic axis ratios of ~ 2/3. Therefore, his 
most prolate models are less elongated than the mean axis ratio predicted by the observations. His prolate 
equilibria are also only obtained when the cores are not very centrally concentrated, with a density contrast 
of only ~ 3 between the centre and the surface. We shall show that our models are both more elongated 
and more centrally concentrated. 

One way that a core could acquire a helical field is by inheriting it from its parent filamentary cloud. 
We previously derived models of filamentary clouds threaded by helical magnetic fields in Fiege and Pudritz 
1999a, (hereafter FP1). We showed that helical fields with constant poloidal and toroidal flux to mass 
loading are consistent with both the virial properties of filaments and the existing data on their density 
profiles (Alves et al. 1998, Lada et al. 1998). Our main goal in this paper is to determine the equilibrium 
states that are accessible to cores that might form from elongated fragments of these magnetized filaments. 
Three parameters are required to specify our filament models, which we discuss in the next section. In 
addition to these three parameters, we only need to specify the mass of a core to calculate its equilibrium 
structure in our model. 

Our main findings are that a large range of core models with projected axis ratios (p) ranging from 
~ 0.3 to ~ 1 can be obtained from our models of filamentary clouds threaded by helical fields. Our 
converged solutions define the prolate magnetic analogues of Bonnor-Ebert (1956, 1955) spheres. Their 
magnetic structure is characterized by a relatively strong "backbone" of poloidal flux surrounded by a 
region where the toroidal field dominates, but usually does not exceed the peak poloidal field. 

We rather completely map out our parameter space by searching for models whose parent filaments 
obey the observational constraints discussed in FP1, and whose fragmentation is driven predominantly by 
self-gravity rather than by MHD sausage instabilities (See Fiege & Pudritz 1999b for details, hereafter 
FP2). The prolate core models that converge are found to occupy a well-defined region of parameter space. 

1 We use the term "helical" rather loosely to describe any field geometry that includes both poloidal and 
toroidal field components. The field lines in our prolate core models are not true helices in the mathematically 
correct sense. 
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The failure of models to converge outside of this region can be traced to non-equilibrium resulting from 
gravitational and MHD instabilities. We also show that the magnetic analogue of the Bonnor-Ebert critical 
mass is somewhat reduced by the helical field in our models. 

A brief outline of our paper is as follows. In Section || we outline the self-consistent field method that 
we use to generate our equilibria. We describe our numerical method in Section ^ and show a gallery 
of representative models in Section ^. In Section |^, we show sequences of models with varying mass and 
external pressure, and show that the Bonnor-Ebert critical mass is reduced by the magnetic field. We derive 
a singular model of isothermal cores wrapped by purely toroidal fields in Section ^, and discuss our results 



Our basic strategy is to use a relatively minor generalization of the self-consistent field method 
originally developed by Mouschovias (1976a), (subsequently used by Mouschovias 1976b, TIN, TIN88b, 
TIN89, TIN90, and Tomisaka 1991) to follow the relaxation of finite fragments of the filamentary clouds 
discussed in FP1. The reader is referred to TIN for a full derivation of the equations. We briefly outline the 
method and discuss our modifications in this section. 

The self-consistent field method iteratively solves Poisson's equation and the equations of 
magnetohydrostatic equilibrium for a self-gravitating cloud threaded by a magnetic field and truncated 
by an external pressure. The source of the external pressure, which is presumably relatively low density 
molecular or atomic gas, is assumed to be non-self-gravitating and of negligible density. Thus, the external 
pressure Ps serves mainly to define the surface of the cloud in our model. 

The self-consistent field method separates the equations of magnetohydrostatic equilibrium into 
equations describing the detailed balance of gravitational, magnetic, and hydrostatic forces parallel to 
the poloidal field lines and perpendicular to them. Following TIN, we use the flux quantity $ which is 
proportional to the true magnetic flux $ mas : 



in Section 0. 



2. 



The Self-Consistent Field Method 



mag 



= 2tt$. 



(1) 



The equation of magnetohydrostatic equilibrium is given by 



Bx (VxB) 

47T 



+ pV* = 0. 



(2) 
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Writing the toroidal field as 

B 4 = ^, (3) 



r 



it follows directly from the azimuthal component of equation [| that is constant along field lines (See TIN 
for proof); hence b$ is a function of $ alone. 

The magnetic field exerts no force in the direction parallel to the poloidal field lines. Thus the equation 
of magnetohydrostatic equilibrium reduces to the hydrostatic equation along a field line: 

— +p— = 0, 4 

OS OS 

where s is the distance along the field line (s = at the symmetry plane defined by z = 0), P is the effective 
pressure, which may contain a contribution due to non-thermal motions of the gas, p is the density, and W 
is the gravitational potential (We do not include the effects of rotation in our analysis.). We assume that 
the gas is "isothermal" in the sense that the total velocity dispersion a is constant: 

P = a 2 p. (5) 

In this case, equation |J can be integrated to obtain 

P-^e-V, (6) 



where Q( < E ) ) is constant along any field line [], but may vary across them. We note that our Q is identical to 
the q used by TIN; we use this different notation to avoid confusing Q with the intrinsic axis ratio q of a 
core. The poloidal components of the magnetic field are related to the flux $ by 

_ 19$ 

Or — 7T- 

r or 

B, . i* (7) 

r oz 

With the help of equations |3| and |6|, the magnetostatic equilibrium equation (equation |^) can be 
re-arranged to give an elliptic type partial differential equation for the azimuthal component of the magnetic 
vector potential (see TIN, equation 2.26): 

where is related to the flux by 

$ 

A* = -, (9) 
r 



The field constant Q is related to the Bernoulli integral along a field line, which is just a 2 InQ. 
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and the differential operator Ai is given by 

A - —+r— (-— ) (10) 
dz 2 dr \r dr J 

We write the left hand side of equation || in terms of A^, rather than <!>, since this is the form of equation ^ 
that we actually finite difference in Section || and solve in Sections || and [|. 

Finally, we require Poisson's equation for the gravitational potential: 

V 2 *=47rGp, (11) 

where p is related to Q and 4" by equation ||. 

It remains only to specify the functions Q(&) and b,p(&) to close the system of equations given above. 
TIN derive the functional form of Q(&) by assuming that cores relax to equilibrium from initially uniform, 
spherical clouds threaded by a constant poloidal magnetic field. We assume instead that cores relax from 
ellipsoidal fragments of a parent filamentary cloud that is threaded by a helical magnetic field (See FP1). 
Therefore, we extend the work of Mouschovias (1976a, b) and TIN by examining the equilibrium states of 
cores that relax from filamentary, rather than spherical, initial clouds. 

The mass contained between poloidal flux surfaces defined by $ and <I> + <5$ is given by 

SM = [ AnrpSw ds, (12) 
Jo 

where s is the distance along the field line, and s m is its maximum value at the surface of the core, where 
the pressure drops to Pg. The thickness 5w of the flux tube varies along its length in such a way that the 
flux <5$ between the flux surfaces remains constant: 

r\B pol \ 

where \B po [\ = B 2 + B 2 is the magnitude of the poloidal magnetic field. Combining this equation with 
equations [| and [l^, we easily obtain the mass to flux ratio: 

where we have used the fact that Q is constant along field lines to take it outside the integral. Solving for 
Q(3>) and replacing the integration over s with an equivalent integral over z, we obtain 

a 2 dM/d^ 



Q(«) 



4tt 



e -BT dz 



(15) 
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where the integral is to be performed along a poloidal field line, and z = z m at the surface of the cloud. 

In order to use equation |l^, we need to evaluate the mass to flux ratio dM/d<& as a function of $ for 
the initial filament, since this function is preserved by flux- freezing as the core settles toward its equilibrium 
state. As discussed in FP1, our models of magnetized filamentary molecular clouds require 3 parameters, 
namely the poloidal and toroidal flux to mass ratios T z and and a third parameter C that determines 
the radial concentration of the cloud. The poloidal and toroidal magnetic fields are given by our filamentary 
cloud model as 

B z = T z p 

= T^rp. (16) 
The surface of the filament, where the pressure drops to that of the ISM, is located at radius 

Rs = r e c . (17) 

where ro is the core radius defined by 

"» = isk (18) 

and p c is the central density of the filament. We consider an ellipsoidal "fluctuation" centred on the axis 
of the filament and assume that all of the gas initially inside of the ellipsoid becomes part of the core in 
its equilibrium state. We assume that the semiminor axis of the ellipsoid is equal to the radius Rs of the 
filament. The semimajor axis Zq (in the direction of the filament axis) is then determined by the mass 
chosen for the core, as we now demonstrate. 

The poloidal field lines threading a filament run parallel to the filament axis since there is no radial 
magnetic field. Therefore, surfaces of constant flux are concentric cylinders centred on the filament axis. 
The mass and flux contained within the ellipsoidal fluctuation between cylindrical radii r and r + Sr are 
given respectively by 



SM = AnrpZo ^l-r 2 /R s 2 Sr (19) 

5$ = T z rp Sr. (20) 

The second of these equations can be integrated to give the total flux $ threading the ellipsoid. We find 
that $ is related to the mass per unit length m of the filament by 

* ci = — . (21) 
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The mass to flux ratio is obtained by combining equations [19] and g0| 

dM 4irZ 



M Tz v l-rV^- (22) 

We note that this equation gives dM/dQ as a function of $ alone (numerically), as required by equation 
|l5| , if we consider r to be a function of $ within the filament. The semimajor axis Zq which appears in 
equation B2| is obtained in terms of the mass and flux by integrating this equation and rearranging: 



Zo = ^kr[J v 1 - 7 - 2 /^ 2 ^, (23) 

which is calculated numerically. 

Ampere's law shows that the function &0($) is directly proportional to the total poloidal current 
passing through the flux surface threaded by flux $. Since b<j, is a function of $ alone, the poloidal currents 
associated with the toroidal field must flow only in the direction of the poloidal field. Tomisaka (1991) 
assumed that tx <J>. However, we obtain the function fo</.( t &) in our model by assuming that cores form 
without dissipation from the initial filamentary state, so that the current flowing between any two flux 
surfaces is a conserved quantity. Specifically, the functional form of bj,($>) follows from equations || and [if]: 

= i> 2 p, (24) 

where r and p are the radius and density of the parent filament, considered to be functions of the flux <E>. 



2.1. Dimensionless Equations 

Following FP1, we write our equations in dimensionless form by normalizing all quantities with respect 
to the velocity dispersion <jfu of the filament and the density p c at its centre. The dimensional scale factors 
can then be written as follows: 

°l = a fa 
Pa = Pc 

I n 



'0 



inGpc 



Pa = Qa = &}uPc 
B a - 
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A = 6^,0 = r B 

$o = rlB . (25) 
The dimensionless forms of the flux to mass ratios T z and are given in FP1 (equation 33) as 

? 2 \ P 

' U ' Ls (20) 



Normalizing all quantities with respect to these dimensional scaling factors, only Poisson's equations 
changes form: 

V 2 * = p. (27) 

The dimensionless forms of the remaining equations are identical to their dimensional forms given by 
equations p[ ^, 0, |^, and [l5|, where (dimensionless) a — 1 is specified as an input parameter ^ , and p, 
P, A^, <I>, B, and Q are all solved self-consistently. 



2.2. Boundary Conditions 

We assume that the poloidal magnetic field threading the original filament is continuous at its surface 
and constant in the external medium, so that A^ oc r in the external medium. The total flux threading our 
grid, of radius r max > R s , is given by 

<t> max = ^ [m + n(r 2 max - R s 2 )P s /a 2 ] . (28) 

Therefore, A <j>tmax is given by 

Affr^rnax ^raaxi (29) 



where c is a constant: 

Gz 



c = 



2?rr 2 

max 



[m + n(rl lax -R s 2 )P s /a 2 ]. (30) 



We find it convenient, for computational reasons, to define a "modified" vector potential A such that 

A = A lf) - cr. (31) 



3 We only consider the case where a is the same in the core as in the parent filament (dimensionless 
(7 = 1). However, we leave a in our equations so that this condition might be relaxed in future work. 



-li- 



lt is easily verified that A satisfies the same equation (equation ||) as A$. However, A vanishes at the outer 
radius of our grid, where r = r max . 

To obtain the boundary condition on A at large z, we assume only that there is no radial component 
of the magnetic field at the edge of our grid corresponding to z = z max The boundary conditions on \P 
and A (in our dimensionless variables) are as follows. 
r = : 

d^/dr = 0, A = (32) 

M , , 

vf= — , A = (33) 

z = 0: 

dV/dz = 0, dA/dz = (34) 



* = dA/dz = (35) 



47T-/ 



3. Numerical Method 

We briefly describe our numerical implementation of the self-consistent field method in this section. 
We work in modified radial coordinates defined by 

r = S r tan £ 

z = S z t&nr) (36) 

so that radial infinity and infinity along the axis of symmetry are respectively mapped to £ = tt/2 and 
i] = tt/2. Our grid is uniformly spaced in £ and rj and contains 2 7 + 1 = 129 points in each dimension. The 
advantage of working in these modified coordinates is that they provide a good dynamical range (spatially) 
for our calculation. We obtain good resolution within the core by keeping a large fraction (typically 
~ 1/3 — 1/2) of the grid cells within its boundary, but our grid still extends to large radius where the 
boundary conditions stated in the previous section apply. 

4 This differs slightly from TIN, who assume that the magnetic field connects to their initially uniform 
field at large z. 
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We begin by initializing our numerical method as follows. We find an equilibrium solution for a 
magnetized filament using the method discussed in FP1. From this solution, we calculate dM/d&, Zq, and 



b,p($) using equations [22|, |23j, and 24 respectively. We begin with a initial guess of the density distribution 
p{r,z) (eg, truncated Gaussian), as well as the magnetic vector potential A(r,z), which gives $(r, z) by 
equations || and |3l|. We then solve our system of equations iteratively as follows: 

1) We finite difference equation |2^, and solve for the gravitational potential 4'(r, z) 

2) Next, we calculate the poloidal magnetic field using equation |?], and accurately locate the surface of the 
core by finding the surface where the pressure drops to Ps- 

3) We integrate along poloidal field lines (level surfaces of 4>) to solve equation ^5] for Q( < f > ). We differentiate 
Q($) to obtain dQ/d§, which is one of the source terms in equation ||. 

4) We interpolate the functions <3($), dQ/d§, and db^/d$ 2 over our grid, using the current value of $ at 
each grid point. 

5) We update A(r, z) by solving the finite difference equation obtained from equation |^, using the functions 
calculated in steps 1 and 5. 

We repeat steps 1-5 until the solution either converges or is rejected. We keep track of the maximum 
relative change in density for all cells whose density is greater than half of the central density. We consider 
a solution to have converged when this quantity either becomes less than 1 x 10 -5 , or it remains below 
1 x 10~ 4 for 50 consecutive iterations. A solution is rejected if the central density falls below that of the 
original filament {p = 1 in our dimcnsionlcss units), or appears to be growing without bound. Our code 
allows a maximum of 750 iterations; this limit is never reached in practice since solutions that converge 
usually do so within ~ 250 iterations. 



4. Exploration of the Parameter Space 

In this section we explore the range of equilibrium models that are accessible to finite segments of 
the filament models described in FP1. We have selected those models from FP1 that agree with the 
observational constraints obtained from our virial analysis of the available observational data (see FP1 
equation 22), and whose kinetic and magnetic energies are nearly in equipartition, with 0.5 < M/K < 2. 
Moreover, we have selected only filaments that fragment slowly, with growth timescales 

/ 77 \ — 

^^ L8 (io^) Myr (37) 
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We refer the reader to FP2 for a discussion of fragmentation in filamentary clouds and the normalization 
used in equation The reason for this selection is to exclude filament models that are dominated by 
rapidly growing MHD-driven sausage instabilities, rather than gravitational instabilities. A total of 71 
filament models from FP1 and FP2 meet these criteria and are considered as initial conditions for our core 
models. 

The Bonnor-Ebert (1956, 1955) critical mass for unmagnetized cores bounded by surface pressure P s is 
given by 

4 

M b .e. = 1.18 — (38) 

Although the ordered magnetic field in our models modifies the critical mass from its hydrostatic value, 
we nevertheless find Mb.e. to be a useful fiducial mass scale to compare with our magnetized models. We 
find that an equilibrium state is likely to be found if the mass of a core is chosen to be comparable to, but 
generally not in excess of the Bonnor-Ebert critical mass for non-magnetic cores. Therefore, we attempt to 
compute core models with masses ranging from 0.3 to 0.8 (in steps of 0.1) times Mb.e. for each of the 71 
filaments chosen as possible initial states. Models that are more massive than ~ 0.8 Mb.e. generally do not 
converge because the central density grows without bound until the numerical calculation is halted. Models 
less massive than ~ 0.3 Mb.e. may converge numerically, but they do not meet our selection criterion that 
the central density of the core should be at least that of the initial filament (See Section ||). Even within 
our chosen mass range, only 152 of the 426 attempted models (~ 1/3) converged and satisfied our condition 
on the central density. 

Figure [j] shows 6 cuts through our parameter space made by holding M/Mb.e. constant. The solid 
dots represent models that converged, while the open circles represent failed or rejected models. It is 
apparent from the figure that most of the converged models have 0.4 < M/Mb.e. < 0.7, with only 4 models 
found with M/Mb.e. — 0.8 and 5 with M/Mb.e. — 0.3. Within each cut of constant M/Mb.e., we find 
that our converged models are separated from those that did not converge by a fairly sharp boundary, with 
little intermixing. This suggests that models that fail to converge likely represent genuine non-equilibrium 
or unstable states. Models do not converge when M ~ 0.8 Mb.e. because the Bonnor-Ebert critical mass 
is reduced in our model by the toroidal field, which works in concert with gravity and surface pressure in 
squeezing the cores. This result is in agreement with Tomisaka (1991) and Habe et al. (1991). We also find 
that r^, never exceeds T z by more than a factor of ~ 1.7 for any converged model. This suggests that MHD 
instabilities might be preventing the convergence of models with very strong toroidal fields, since we found 
in FP2 that MHD sausage instabilities were triggered in filaments when T$/Y z ~ 2 (See also discussion in 
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McKee 1999). 

4.1. A Gallery of Models 

We have computed a total of 152 models scattered throughout our parameter space. Figure || shows 
a gallery of representative models, with a range of intrinsic axis ratios q ranging from 0.2 to 0.8, shown in 
order of increasing mass. The density contours are linearly spaced from the central density (given in each 
panel) to the density at the surface of pressure truncation. We find that the density contours near the centre 
of our models are nearly circular, but become increasingly elongated towards the surface. This is explained 
by the dominance of gravity near the centre, and the increasing importance of the toroidal magnetic field 
in the outer regions. The density contours in our models are therefore somewhat more complex than in the 
empirical model used by Myers et al. (1991) and Ryden (1996), who assume that the density contours of 
cores can be described by nested ellipses of constant axis ratio. 

Our models have a central "core" region of nearly uniform density, surrounded by an outer "envelope," 
where the density gradually decreases with radius. Figure ^ shows that our models are generally quite 
truncated by the external pressure, typically extending to only ~ 2ro in cylindrical radius. Therefore, they 
are truncated before they attain the r~ 2 density profile (or any other power law) that characterizes the 
outer envelope of a Bonnor-Ebert sphere. 

For each of our models, we have computed the intrinsic axis ratio q by calculating the best-fitting 
ellipse (in a least squares sense) to the half maximum density contour. We have also calculated the mean 
projected axis ratio (p) for each model using the procedure outlined in Appendix A. (Essentially, (p) is 
the expectation value of the projected axis ratio of the half maximum surface density contour for a core 
oriented at a random inclination angle relative to the plane of the sky.) 

We show histograms of the intrinsic and mean projected axis ratios in Figure ^ below, where we find 
that 0.1 ~ q ~ 0.9 and 0.3 ~ (p) ~ 1. These histo grams should not be interpreted as a prediction of the 
distribution of axis ratios for a sample of cores, as shown in Ryden (1996) for example. What they do 
show is the distribution of axis ratios that was actually obtained by our sample of converged models, which 
provides an indication of the relative ease or difficulty in producing core models with a given axis ratio. 
The mean value of q derived from observations is probably between 0.4 and 0.5 (Myers et al. 1991, Ryden 
1996), but the predicted distributions are quite broad, containing significant numbers of cores that are 
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much more elongated or more nearly spherical (Ryden 1996). Therefore, Figure^ shows that our models 
cover the likely range of axis ratios. 

Figure ^ shows that the relationship between (p) and q, numerically determined by our model, is in 
reasonable agreement with the analytical relation (dashed line) used by Myers et al. (1991, their equation 
1), even though the internal structure of our models is quite different from theirs. The relation between q 
and (p) is clearly not very sensitive to the detailed internal structure of cores. The error bars represent the 
standard deviation in (p) for a sample of cores with axis ratio q, oriented randomly relative to the observer. 

We also estimate the degree to which the half maximum density and surface density contours of our 
models deviate from a pure ellipse. Our method for this analysis is the Fourier method developed by 
Bender and Mollenhoff (1987), who used the method to quantify the "diskiness" or "boxiness" of elliptical 
galaxy isophotes. Their method works by finding the best-fitting (in a least squares sense) ellipse to an 
isophote and Fourier analysing the residual. Their analysis shows that the fourth Fourier component 
of the residual dominates the Fourier spectrum, and when normalized with respect to the semimajor axis 
a provides an excellent indication of the deviation from an ellipse. For our purposes, a^/a > describes 
density or surface density contours that are more peaked than an ellipse along the axis of symmetry (ie: 
shaped like a North- American football or lemon), while a\ja < describes contours that are more flattened 
than an ellipse. 

For each of our models, we calculate both an intrinsic a^/a parameter for the half maximum density 
contour, and a mean projected {a^/a) parameter for the half maximum surface density contour, averaging 
over all possible orientations using the procedure described in Appendix We show the distributions 
of these shape parameters for our models in Figure |. We find that all of our models have intrinsic and 
projected shape parameters that are greater than or equal to zero, but usually less than ~ 0.05. Therefore, 
many of our models are slightly football-shaped. This might provide one way to distinguish observationally 
between our models and flattened models threaded by purely poloidal fields (eg. TIN). 

We find no clear correlation between 04/0 and q. However, very few elongated models (with (p) ~ 0.6) 
are significantly peaked {(a^/a) ~ 0.02). Conversely, those with a high value of a^/a are never very 
elongated. 

An interesting feature found in many of our models is that the field lines bulge slightly outward around 
the core (See Figure || for several examples.). This is opposite the hourglass shaped field lines found by 
Mouschovias (1976b), TIN, TIN88b, TIN89, TIN90, and Tomisaka (1991). Physically, this field structure 
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is mainly due to the radial pinch of the toroidal field, which most significantly affects the density and field 
line structure near the tips of cores (where the symmetry axis intersects the surface of the core) . The gas 
pressure is lowest and the gravitational forces are weakest near the tips, so the toroidal field can squeeze in 
further than in the midplane. 

We find a continuum of models that range considerably in shape and ellipticity. We discuss two 
representative models chosen because of their contrasting behaviour. The model shown in Figure ^| is 
dominated by the B z component of the magnetic field nearly everywhere. The toroidal field vanishes at the 
axis of symmetry and increases radially outward, becoming comparable to the poloidal field only very near 
the edge of the core. It is strongest near the tips of the cores rather than in the midplane. The radial field 
component is much smaller than B z everywhere within the core. Many models like the one shown in Figure 
||, with a strong "backbone" of poloidal flux, are very elongated. We find models similar to the one shown, 
but with intrinsic axis ratios as low as ~ 0.1 and mean projected axis ratios as low as ~ 0.3. The density 
contours are usually quite nearly elliptical when a strong poloidal field is present, with < a^/a ~ 0.03 in 
most cases. 

The model shown in Figure |^ contains regions where the toroidal field becomes significantly stronger 
than the poloidal field. The density contours are quite peaked towards the axis of symmetry, with 
di/a > 0.04 for some models like the one shown. While some models with this sort of behaviour are fairly 
elongated, it is interesting that most are actually less elongated than those with strong poloidal fields. 

4.2. Trends in the Model Parameters 

Figures ^ and |^ illustrate that our models generally become more elongated as the mass is decreased or 
as Ps/(P) is increased. The reasons for these trends are straightforward. Given that rotation is negligible 
and the velocity dispersion is isotropic, gravitationally dominated cores should be spherical. We find nearly 
spherical, gravitationally dominated cores when the mass is large and the central density is high. The 
helical magnetic field can most effectively shape low mass cores, in which self-gravity is less important and 
the surface pressure is a significant fraction of the central pressure. The most massive core in our sample 
of converged solutions is shown in Figure |l(j. We show the radial density structure of a purely hydrostatic 
isothermal sphere overlaid on radial and axial cuts of the density from our massive solution. We find that 
the inner, most gravitationally dominated regions are described by an isothermal sphere to a high degree 
of accuracy within ~ 0.5 of the centre (maximum deviation = 1.3%). Beyond this radius, the magnetic 
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field becomes important relative to gravity, resulting in an axis ratio of 0.64 for the contour describing the 
surface of the core. 



Figure 11 verifies our assertion in Section 4.1 that only modest toroidal fields are required to produce 
substantially elongated cores. In panel (a), we show that the projected axis ratio (p) decreases as the ratio 
Mtot/K increases, where M tot and K are respectively the total magnetic energy ^] and the kinetic energy, 
defined by 

Mtot = ^- [ B 2 dV 

K ee lJ y PdV. (39) 

Models with (p) < 0.6 have Mtot/K ~ 0.4 in our model, but Mtot/K <~ 2 even for our most elongated 
models. Therefore, the field strengths required by our model are consistent with approximate equipartition 
between kinetic and magnetic energies, which has often been inferred from observations (eg. Myers & 
Goodman 1988a, b). Panel (b) shows that most of the magnetic energy is in the poloidal component of the 
field for significantly elongated cores, with (p) < 0.6. We plot (p) against M^/M po i, where M^ and M po i 
are respectively the toroidal and poloidal contributions to the total magnetic energy in equation g9[ As a 
general trend, we find that the most elongated models have the lowest values of M^/Mpoi, and that most 
models with (p) < 0.6 have 0.3 ~ M$/M vo i ~ 1. This trend might seem paradoxical because the toroidal 
field is responsible for the elongation. However, models with energetically significant toroidal fields must 
also have strong poloidal fields to suppress MHD instabilities. 



5. Sequences of Models with Varying Mass and External Pressure 



Our emphasis in Sections 4T and 4^2 was on exploring the full range of equilibrium states that are 
accessible to our models. In this section, we select a few filament models and examine sequences of core 
models that are obtained by gradually increasing the mass and the external pressure. 

The sequences shown in this section terminate when the mass or pressure become sufficiently high that 



5 The magnetic energy M to t is not the full magnetic energy term from the virial theorem, since we do not 
include any surface terms. Our definition is more relevant to Zeeman measurements, since M tot is related to 
the total average field strength within the core. We also note that Mtot/K is analogous to our parameter X 
in FP1. 
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no equilibrium can be found due to a central density that increases with each iteration of our numerical 
code, apparently without bound. The maximum mass or pressure found in this way clearly represents the 
Bonnor-Ebert critical mass or pressure, which are both reduced by the magnetic field. These findings are 
consistent with the results of Tomisaka (1991) and Habe et al (1991). 



5.1. Effect of Mass 

Figure |l2| shows how the intrinsic axis ratio q varies with mass for core models obtained from 3 different 
parent filaments. We have normalized the mass by the Bonnor-Ebert critical mass Mb.e., to provide a 
convenient reference point for the mass. The filament parameters T z , F^, and C are given in Table [l], 
as well as the external pressure P ex t and Mb.e.- We find that the axis ratio q decreases with increasing 
mass, when the mass is very low (less than a few tenths of Mb.e.)- These weakly self-gravitating equilibria 
are strongly pinched by the toroidal field and more closely resemble short filaments than cores. Once 
self-gravity becomes important (when the mass is a significant fraction of Mb.e.), the axis ratio decreases 
with increasing mass. Our models become nearly spherical when M core ~ 0.75A1b.e.- 

The sequences shown in Figure |l2| all terminate at masses that are about 20% less than Mb.e. , although 
we have attempted to compute models with greater masses. This is consistent with our finding that very 
few models converge when M cr i t /MB.E. > 0.8 (see Figure [TJ). Therefore, the critical mass M cr u < Mb.e. 
for our models, as a result of the pinch provided by the helical magnetic field. 



Sequence 


r 2 




C 


Pext 


Mb.e. 


1 


7.17 


6.05 


0.165 


0.436 


79.6 


2 


6.76 


6.45 


0.200 


0.353 


88.6 


3 


4.53 


6.18 


0.235 


0.243 


106.6 



Table 1: We give the dimensionless parameters used to compute the sequences of models shown in Figure 
|l2| . Parameters to the left of the double vertical line are input parameters, while the quantities to the right 
are derived parameters of the model. 
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5.2. Bonnor-Ebert Stability and the Effect of the External Pressure 

We have previously considered the external pressure bounding the core to be the same as that bounding 
the parent filament. However, for these sequences only, we vary the external pressure to determine how 
squeezing prolate cores affects their shapes. This procedure also allows us to determine the Bonnor-Ebert 
stability of our models. 

Figure |l3| shows how the shape of a core with fixed mass, T z , Tj,, and C changes as the external 
pressure is varied. We have held the mass fixed and varied the external pressure, which we have normalized 
by the Bonnor-Ebert critical pressure Pb.e., defined by 

Pb.e. = 1-39^5- (40) 

core 

for convenience. The parameters describing the parent filament are given in Table |[ as well as the mass 
M core and Pb.e.- We find that our models are most elongated when the surface pressure Ps is well below 
the Bonnor-Ebert critical pressure Pb.e.- Cores become strongly self-gravitating, and therefore nearly 
spherical, when Ps is comparable to Pb.e.- 

We find that all of the sequences terminate at pressures that are well below the Bonnor-Ebert critical 
pressure. Therefore, the critical pressure P cr it < Pb.e., as a result of the helical field squeezing our models 
in concert with the external pressure. 



6. A Singular Model for Prolate Cores 

We have found a very simple singular isothermal solution for a prolate core containing a purely toroidal 
magnetic field. While a purely toroidal field is admittedly not very realistic, the solution is interesting 
because it probably represents the simplest possible model of a magnetized prolate core, and is in fact a 



Sequence 


T z 




C 


M core 


Pb.e. 


1 


7.17 


6.05 


0.165 


39.8 


0.519 


2 


4.53 


6.18 


0.235 


42.6 


0.496 


3 


6.76 


6.45 


0.200 


44.3 


0.484 



Table 2: We give the dimensionless parameters used to compute the sequences of models shown in Figure 
|l3| . Parameters to the left of the double vertical line are input parameters, while the quantity to the right is 
a derived parameter of the model. 
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generalization of the well-known singular isothermal sphere. From a physical standpoint, this model is 
useful because it demonstrates that a relatively moderate toroidal field, with plasma /3 > 1 everywhere, can 
result in elongated cores with axis ratios that are in good agreement with the observations. 

We assume a logarithmic gravitational potential, in cylindrical coordinates (r, <fr, z), of the form 

^ = ^ln (r2+ 2 cV) (41) 

where e is the axis ratio of the isopotential surfaces, and a is assumed to be constant []. This is the potential 

of a singular isothermal sphere when e = 1, and describes singular prolate models when e < 1. In Appendix 

[b|, we show that the corresponding density and toroidal magnetic field are given by 

o 2 e 2 [r 2 + (2-e 2 )z 2 ] 
' 2wG (r 2 + e 2 z 2 ) 2 1 > 



2a 2 rVT^ 

- G l/2 ( r 2 +e 2 z 2)' ^ 

where all quantities are the same as in our numerical models (See Section ^|) . It is clear that more elongated 
solutions (smaller values of e) are consistent with larger values of B^. 

We note that choosing e = 1 recovers the well-known singular isothermal sphere solution, for which B^ 
vanishes. It is interesting that the density vanishes when e = 0. Therefore, the solution does not reduce to a 
singular model for filaments in this limit (See also FP1). An important feature of equation ^ is that there 
is no solution when e > 1. Hence, there are no oblate solutions that are consistent with our singular model. 



The plasma f3 parameter is given by 

(2-e 2 )z 2 



(44) 



which has a maximum of e 2 /(l — e 2 ) in the midplane, where z = 0. Therefore, e = \j\J~2 results in a 
moderately magnetized model with (3 >1 everywhere, which we show in Figure [l4]. We find that the density 
contours of this singular model have an axis ratio of ~ 0.46 (See Section [I] and Appendix |A| for a discussion 
of our method of estimating axis ratios.). This result is in good agreement with the mean intrinsic axis 
ratios that are typically derived for samples of cores (See Ryden 1996). 

A simple numerical example illustrates that relatively moderate toroidal fields can result in significant 
elongation. We consider the magnetic field at a distance of 0.1 pc from the centre of the e = l/y/2 model 



6 See Binney and Tremaine (1994) equation 2-54a for a similar potential used in the context of galactic 
disks. 
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discussed above. This distance corresponds to the radius at which the thermal and non-thermal velocity 
dispersions are approximately equal for low mass cores (Fuller & Myers 1992, Casseli & Myers 1996). 



Therefore, a ~ 0.26 km s , assuming a temperature of 10 K for the gas. Equation 43 requires a maximum 
strength of only s» 12 /iG for at this distance. 

The toroidal field in the singular model is strongest in the midplane, which squeezes the gas so that the 
constant density surfaces are slightly pinched in. This is markedly different from our non-singular numerical 
solutions, which have their maximum toroidal field off the midplane. However, models with purely toroidal 
fields, such as our singular solution, are a separate class of solutions that cannot be obtained from the 
self-consistent field equations (discussed in Section ||) , even in the limit of a vanishing poloidal field. The 
reason for this difference is that azimuthal component of the magnetohydrostatic equilibrium equation 
(equation ^) is automatically satisfied for purely toroidal field geometries. In the more general case of a 
helical field, torsional equilibrium requires 6^ to be constant along poloidal field lines, so that the class of 
allowed solutions is somewhat more restricted. 



7. Discussion and Conclusions 

We have previously (FP1) suggested that some filamentary molecular clouds may contain helical 
magnetic fields, and that such filaments are unstable against fragmentation along their lengths (FP2). Some 
fragments may dynamically evolve into equilibrium states that can be identified as cloud cores. The central 
point of the present work is to demonstrate that most of the possible equilibrium states that arise from our 
filament models are prolate, and that many have realistic axis ratios. 

It is very encouraging that we are able to explain several of the observed features of filamentary 
clouds and their cores within a self-consistent theoretical framework. Future tests of our models may be 
provided by measurements of polarized dust grain emission. A useful test of our models will be to predict 
the polarization due to both our filament models (FP1) and the prolate core models derived in this paper. 
This is done in a forthcoming study. Perhaps even more importantly, the central assumption of our model, 
which posits a constant flux to mass loading, requires direct confirmation. This could one day be feasible 
with sufficiently high spatial resolution Zeeman measurements. 

Previous observations of the shapes of cores have concentrated mainly on their axis ratios (Myers et 
al. 1991, Ryden 1996). However, we have found that the shapes of cores can usefully be described by a 
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second observable parameter 04/0, which quantifies the deviation of the density contours from ellipses. The 
toroidal field in our model results in a^/a > 0, which is consistent with surface density contours that are 
somewhat football-shaped (ie. peaked along the major axis). Therefore a detailed isophotal analysis of 
cores may provide important clues to their internal magnetic structures. 

Core models that exceed the critical mass for stability in the sense of Bonnor (1956) and Ebert (1955) 
are doomed to collapse. The dynamical collapse of one of our prolate core models from a critical initial 
state is a difficult problem that should be addressed by future numerical simulations. 

8. Summary 

1 . Most of the core models that are consistent with our previous models of filamentary molecular clouds are 
found to be prolate. A wide range of axis ratios is possible, with 0.1 ~ q ~ 0.9 and 0.3 ~ (p) ~ 1. 

2. Both density and surface density profiles are found to be slightly football-shaped (ie. peaked toward the 
symmetry axis) in our models, as quantified by the 04 /a parameter of Bender & Mollenhoff (1987). 

3. The Bonnor-Ebert critical mass is reduced by about 20% in our models by the toroidal field. 

4. We find that models are generally more elongated for low values of the mass ratio M/Mb.e. and the 
pressure ratio Ps/Pb.e.- 

5. Only modest toroidal fields are required to produce elongated cores. Most of our models have < B z 
nearly everywhere. 
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A. The Mean Projected Axis Ratio and a^ja Shape Parameter 

For each of the models discussed in the previous sections, we have calculated the surface density at 
various angles of inclination. By assuming that cores are oriented randomly, we calculate a distribution 
function for the axis ratio and the a±/a shape parameter of the half maximum surface density contour, from 
which we obtain their mean values and standard deviations. We work in Cartesian coordinates (x 1 , y', z'), 
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with the origin at the centre of the core, the y' axis along the line of sight to a distant observer, and the 
x' — z' plane parallel to the plane of the sky. We define a second set of Cartesian coordinates (x, y, z) fixed 
relative to the core, with z corresponding to the axis of symmetry. The z axis is assumed to be tilted by 
an inclination angle i relative to the plane of the sky and towards the observer. The coordinates are then 
related by 



x = x 



y — y cos i — z sin i 

z' = ysini + zcosi. (Al) 



Since the cylindrical radius (relative to the symmetry axis) is given simply by r = \J x 2 + y 2 , we interpolate 
p{r,z) over the primed coordinates to obtain p(x',y',z'). The surface density Yl(x',z') is then obtained 
by integrating over y' . We then compute the contour of the half maximum surface density, from which 
we estimate the axis ratio by fitting an ellipse as described in Section ^. We also compute the a^/a shape 
parameter, using the method discussed in Bender and Mollenhoff (1987). We repeat this procedure at many 
different inclination angles to determine the projected axis ratio p and a^/a as a function of i. 

The z axis is assumed to be oriented in a random direction (#, <p) relative to the z' axis. The inclination 
angle is related to 9 and <j> by 

cos i = 1 - sin 2 9 sin 2 4>. (A2) 
The probability dp, of the z axis being oriented between <f> and <f> + d<j) is 

dp, = (A3) 

Integrating, the probability of the z axis being oriented between and <fi is given by </>/27r, so that 

= 2tt£ (A4) 

if £ is a uniformly distributed random variable. Similarly, it is easy to show that 

cos6l = l-2£ (A5) 

for another random variable £. We generate a large number of realizations (typically 10 5 ) of the orientation 



from which we calculate the orientation i by equation A2. We interpolate to find p and 0,4/ a for each 



realization, and calculate the mean values and standard deviations of the resulting distributions. 
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B. Derivation of the Singular Isothermal Model for Prolate Cores 

Working in the dimensionless variables introduced in equation |2^ and assuming that the magnetic field 
is purely toroidal, the equation of magnetohydrostatic equilibrium (equation [2]) can be written as 

87rr 2 e~*VQ = \Ib\ (Bl) 

where b^ is related to B,p by equation The variable Q has the same form as in equation ^, but is no 
longer contsrained to be a function of 

Q = pe s . (B2) 



We note that the r and z components of equation Bl can be combined into the form of a Jacobian 



determinant, which implies that 6^ can be considered as a function of Q alone. 

The potential for the singular isothermal sphere model, in our dimensionless variables, is given by 
^ = ln(i? 2 /2), where R is the spherical radius. Therefore, we consider distorted logarithmic potentials of 
the form 

r 2 + e 2 z 2 



vf = In ^ j (B3) 

where e is the axis ratio of the isopotential surfaces. The density is easily obtained from Poisson's equation 
(equation |7|) : 

Jr 2 + (2-e 2 )z 2 l , 
P = 2e 2 1 2 I 2 ' J , B4 

which appears in dimensional form as equation Wa. 



Inserting equations B2 and B4 into equation |B2], wc find that Q can be written in the form 

£ 2 [l + (2- £ 2 )C 2 ] 
1 + e 2 C 



Q(0 = 1 ; , J , (B5) 



where 

C = z/r. (B6) 



Since is a function of Q alone, we can consider both Q and b$ to be functions of £ in equation Bl. This 
equation becomes an ordinary differential equation in for b$(£): 

db^ 8?r dQ 

dC l + e 2 C 2 d(' 1 ' 

Substituting Q(C) into this equation and integrating, we find that 



2 16tt(1 - e 2 v 

h= ( l + e 2 C 2)2- (B8) 
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The toroidal field is directly obtained from this equation by replacing both £ and 6^ with their respective 
definitions (equations B6 and|^). The final dimensional form is written out in equation |43|. 
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M/M D £ =0.4 





Fig. 1. — We show the regions of parameter space where our models converge (black dots) or fail to converge 
(open circles). Each attempted solution is consistent with a parent filament obeying the observational 
constraints established in FP1 and having rather low fragmentation growth rates, as given by equation |37|. 
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r z = 6.46 
r = 8.23 

C = 0.112 
M/M B E = 0.4 

M = 36.4 
P s = 0.333 

P = 113 
q = 0.135 
<p> = 0.294 
a Ja = 0.0272 



1 r 



4 

r/r„ 



r z = 8.06 
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46.5 
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q = 0.285 
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4 
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C = 0.165 
M/M c = 0.6 



M = 47.8 
P s = 0.436 
P c = 2.75 

q = 0.535 
<p> = 0.664 
a 7a = 0.00649 



4 




Fig. 2.- 
the range 



We show density contours and poloidal field lines for a set of representative models to illustrate 
of equilibria that we obtain. Models are shown in order of increasing q. 
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Fig. 3. — We show the distribution of intrinsic and mean projected axis ratios obtained for our converged 
models. All 152 models are elongated, many with axis ratios in good agreement with the observed values. 
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Fig. 4. — We show the relationship between the intrinsic and mean projected axis ratios for our model. 
We find reasonable agreement with the formula (dashed line) given by Myers et al. (1991) for an empirical 
model of prolate cores. The error bars represent the standard deviation of p about (p), assuming that cores 
are oriented randomly on the sky. 
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a 4 /a (intrinsic) 




0.05 0.1 
<a /a> (projected) 



Fig. 5. — We show the distribution of the intrinsic and mean projected cla/cl shape parameters for our 
converged models. 
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Fig. 6. — We show the structure of a representative model whose toroidal field is relatively weak compared 
to its poloidal field nearly everywhere, panel (a) (upper left): density contours and poloidal field lines, 
panel (b) (upper right): 

Contours of toroidal field intensity, panel (c) (lower left): Density cuts taken radially in the midplane (solid 
line) and along the asymmetry axis (dashed line), panel (d) (lower right): The variation of B z (solid line), 
Bff, (dashed line), and B r (dot-dashed line) taken along a radial cut in the midplane. 
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0.6 

M/M B.E. 




0.6 

M/M B.E. 



Fig. 8. — We show how the intrinsic and mean projected axis ratios depend on the ratio of the mass to the 
critical mass Mb.e. for non-magnetic cores. Gravitationally dominated cores with masses near Mb.e. are 
much less elongated than those with smaller masses. The error bars on p are suppressed for clarity. 
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Fig. 9. — We show how the intrinsic and mean projected axis ratios depend on the ratio of the external 
pressure to the pressure at the centre of the core. Pressure dominated cores with relatively large values of 
Ps I {P} are more elongated than less pressure dominated models. 
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Fig. 10. — We show the most massive core obtained in our sample ol 152 converged solutions. The heavy solid 
line and dashed line represent density cuts taken radially in the midplane and along the axis of symmetry 
respectively. Both cuts agree well with an isothermal sphere (thin solid line) in the inner, gravitationally 
dominated region, but deviate significantly in the pressure dominated region near the surface. 
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Fig. 11. — Panel (a) (left): we show how the projected axis ratio (p) depends on the ratio of the magnetic 
to kinetic energy M tot /K. M to t does not exceed <~ 2K for any of our models. Panel (b) (right): We 
show how (p) depends on ratio of the toroidal to poloidal magnetic energy M<p/M po i. Significantly elongated 
models have < M po i. 
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Fig. 12. — We show the effect of varying the masses of cores on q while holding all other parameters fixed. 
The curves are labeled to correspond with the column labeled "Sequence" in Table 111. 
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Fig. 13. — We show the effect of varying the external pressure on q for cores of constant mass. The curves 
are labeled to correspond with the column labeled "Sequence" in Table @. 
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Fig. 14. — We show density contours (top, panel a) and toroidal field intensity contour (bottom, panel b) 
for the singular model for prolate cores with a purely toroidal field. 
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